Expectation Values of Local Fields in an Integrable Theory after a Quantum Quench 
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The expectation values of local fields of any interacting quantum theory after a quench process 
are key quantities for matching theoretical and experimental results. For quantum integrable field 
theories, we argue that they can be obtained by a generalization of the Leclair-Mussardo formula 
and a Bethe Ansatz result of Caux and Konik. Specializing to the Sinh-Gordon model and taking 
the non-relativistic limit, one can recover the results of Kormos et al. for the Lieb-Liniger model. 
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The aim of this paper is to set up a formula that al- 
lows the computation of the asymptotic time Expectation 
Value (EV) of local fields of a (1 + 1) dimensional Quan- 
tum Integrable Field Theory (QIFT) after a quantum 
quench, i.e. an abrupt change of the parameters of the 
Hamiltonian. The subject of quantum quenches have re- 
cently attracted a lot of attention, both from experimen- 
tal and theoretical point of view, see for instance [1-18] . 
QIFT's are special quantum many-body systems since 
they have an infinite number of conservation laws and 
their dynamics is strongly constrained (see, e.g. [19] and 
references therein): their asymptotic time regime is ex- 
pected to violate the ergodicity property although it may 
be described by a Generalized Gibbs Ensemble (GGE), as 
advocated in [14] and shown in a series of papers, among 
which [15-17]. The many questions coming from this gen- 
eral context (e.g. the setting of equilibration without a 
proper thermalization[l], the importance of rare states 
[12, 13], the phenomena of pre-thermalization [8-11], the 
role of the local as well as the non-local charges and the 
identification of proper density matrix [16], the relation- 
ship between pre-quench and post-quench operators [18], 
the hints coming from Bethe Ansatz [15, 26, 27], etc.) 
have stimulated an intense research activity and a lot of 
progress has been made on many of these topics. How- 
ever, despite these important advances, it seems that a 
proper formula for computing these asymptotic EV in 
QIFT is still missing in the literature. This formula is 
proposed here and illustrated with the avail of the Lieb- 
Liniger model, i.e. the non-relativistic limit of the Sinh- 
Gordon (ShG) theory [23]. It is also worth to mention 
that, if one restores a h dependence, the limit h — > of 
the formula given below identifies the EV in the purely 
classical dynamics of relativistic integrable models [25]. 

In a nutshell, the formula is just a generalization of 
the LeClair-Mussardo formula previously obtained in the 
context of QIFT at a finite temperature [20]. Namely, 
let 0(x,t) be the local field and (?/'o|C(OI 1 /'o) its expec- 
tation value on a macroscopic state \ipo) of the theory 
(i.e. made of an infinite superposition of multi-particle 
states), which we assume to be translation invariant and 
not an eigenstate of the Hamiltonian. For simplicity we 
will consider a QIFT made of only one type of particle, 



the generalization to more general case being straightfor- 
ward. Let's define the Dynamical Average (DA) of this 
field on \ipo) as the time average after the quench at t = 



(O) 
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lim 



dt (ih\o{t)\i>o) 
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Let's instead define the Ensemble Average (EA) of this 
observable on this state as 
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where f{9 l ) = {e^ e ^ -S'(O))- 1 (with S(9) the exact two- 
body S-matrix of the model), while \ 9) = |#i, . . . ,9 n ) 
((0\ = (9„, . . . , 9i\) denotes the asymptotic states of 
the integrable theory expressed in terms of the ra- 
pidities di, with E{0) = m cosh 9, p(9) = msinh.9. 
This formula employs both the pseudo-energy e(6), so- 
lution of a Generalized Bethe Ansatz (GBA) equation 
based on the state \ip ) [27], and the connected di- 
agonal Form Factor (FF) of the operator O, defined 

as (t\0\7) conu = FP (lirn^o W>\t , t - m + ztf>) 
where rf = r) n , . . . ,rji and FP in front of the expression 
means taking its finite part, i.e. omitting all the terms 
of the form r/i/rjj and l/ryf where p is a positive integer. 
Our claim is that it holds the following equality between 
the two averages 
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This equality provides an efficient way of computing the 
EV in a QIFT in terms of a generically very fast con- 
vergent series. As shown below, the Ensemble Average of 
eq.(2) is nothing else but the formula that allows us to 
extract the regularized finite expression of the EV com- 
ing from the Diagonal Ensemble, reached by using the 
Dynamical Average of eq.(l). 

To spell out the content of (3), let's start by showing 
the computation of the EV of local fields in the simplest 
QIFT, i.e. the free theory. In the infinite volume, the 



solution of the cq. of motion (□ + to 2 ) <j>(x, t) = is 
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The modes A{9) and A^(9) are fixed in terms of the 
boundary conditions at t = 0, i.e. by the expectation 
value of the field and its time derivative on the states \ijjo), 
and satisfy the commutation relation [A{9) , A^ {9')] = 
2tt5(8 — 9'). The dynamics of the model is supported 
by the infinite number of non-local conserved quantities 
given by the mode number occupations N(6) = =L \A(9) | 2 
for any 9. To find the infinite set of local conserved quan- 
tities it is convenient to go in the light-cone coordinates 
t = t + x and a = x — t, where the eq. of motion 



becomes 



and, as a consequence, there is 



the infinite chain of conservation laws (n 
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where n(T = d"(f> and analogously for 4> nT . These equa 
tions are of the general form d T A = d a B and, going 
back to the coordinates (x,t), they become the conti- 
nuity equation d t (A + B) = d x (B — A). One can easily 
identify the even and odd charges £271-1 an d O-m-x (the 
first ones are the energy and momentum of the field): up 
to normalization, they can be written as 



£2n-l = "I 
02n-l = 
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They imply that each particle state \9) (as well as all 
multi-particle states) is a common eigenvectors of all 
these conserved quantities, with eigenvalues 



Sn 



to 2 ™- 1 cosh[(2n - 1)0] \9) 



O n \9) = to 2 



sinh[(2?i - 1) 



(5) 



It is pretty evident that the knowledge of the mode occu- 
pation |A((9)| 2 fixes all the local charges but, under both 
general mathematical and physical conditions, it is also 
true the viceversa, alias the knowledge of £2^-1 's and 
£>2n-i's is enough to fix \A{9) | 2 . Being linearly related 
one to the other, the two types of conservation laws are 
then essentially interchangeable. It must be stressed that 
eq. (5) holds also in interacting QIFT as the ShG model, 
for instance, where the only substitution is |^4(#)| 2 with 
the proper action variable P(9) in (5). 

In the free theory, the exact solution (4) of the eq. of 
motion allows us to easily compute the stationary EV 
of any local function F[cj>(x,t)], defined with a proper 
normal-ordering of the operators. Since we are interested 
in field configurations with finite energy density, our the- 
ory has to be defined on a circle of length L and then 
send L — > 00 so that E/L is always finite, even in this 
limit. The momenta of the particles will be quantized 
in unit of 2n/L which become dense when L — > 00. The 




FIG. 1: On the left: time evolution of (ip \ : <j> 2 (t) : \ip ) for 
t = [0, 10 4 £o] and a typical translation invariant state |V>o)- 
On the right: its time average. 



Dynamical Average of cq.(l) smooths out all microscopic 
fluctuations of the obscrvables and leads, asymptotically, 
to stationary values. One way to visualize these EV is to 
consider a quench in the free theory where the mode oc- 
cupations are very large, so that the theory can be treated 
classically and solved in this way [24, 25]: the typical time 
evolution of a space- average observable F[(j>(t)] is shown 
on the left plot of Figure (1) where, as time goes by, there 
are persistent fluctuations. These can be flatten by tak- 
ing the time average: the mean then rapidly converges to 
the asymptotic value (not necessarily a thermal one). 

Let's consider a quench protocol whose instances have 
in common just the same energy density E/L. Since the 
energy is a very degenerate observable, each quench pro- 
cess corresponds to different values of all other conserved 
charges and the typical result is shown in Figure 2, where 
one can see a large spread of these EV. This strong de- 
pendence on the initial data is easily explained. Focus, 
for simplicity, the attention on the EV of the observables 
: (f) k (x,t) : of a free theory. While (: 4> 2n+l :) DA vanish 
by symmetry, a phase-stationary argument for the even 
ones leads to 



(: 4> 2 -)da = 



IDA 
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= (2n-l)!!b n (7) 



Since they explicitly depend on the initial condition 
through the |^4(#)| 2 's, they can never be derived by a 
Gibbs Ensemble average 
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) DA ± Z- 1 Tr(: 
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Notice, however, that the VE (6) and (7) have pre- 
cisely the form of the Ensemble Average, since the non- 




FIG. 2: Time average of (tpo\ : <f> : (t)\ipo) for different initial 
states with the same value of the energy density (in this case 
E/L = 100). 



zero FF arc (0j : 4> 2n (0) : |^) con „ - 2 n (2n)\5 n , n 
while the pseudo-energy is given by e(0) = log(l + 
\A(f3)\~ 2 ), easily derivable for the free theory (with an 
S- matrix S(9) = 1) by using the GGE in the form 
PGGE = Z- l eiq>[-J^e(6)\A(e)\ 2 }. The lagrangian 
multipliers e(0) (one for each mode) are fixed by the 
initial occupation numbers as \A{6) | 2 = (e c ^ — 1) _1 
[5]. The statistical weight given by the GGE is equiva- 
lent to {A(9)tf{9')) EA = 2ttS{9 - 0'), {A{9)A{9')) EA = 
(A* (9)A^ (9')) ea = , from which one can easily get the 
generating function 



(exp[ia</)]) 



EA 
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Expanding in power series in a the left/right hand sides 
and comparing equal powers in a, one recovers the pre- 
vious results (6), (7) . Notice that, while the EV of var- 
ious observables explicitly depend on the initial condi- 
tions, for the free theory there are however some ra- 
tios completely independent from them as, for instance, 



Finite Part of the Form Factors (which are regular func- 
tions of the rapidities) and the pseudo-rapidity e(0) (so- 
lution of a Bcthc Ansatz equation and regular function 
as well) . All these steps are implemented by defining the 
theory on a finite length L and then finally send L — > oo. 
In summary, the finite and regularized expression that 
has to be associated to the Diagonal Ensemble is the one 
given in eq.(2). 

For the quench situation in a QIFT, the only ques- 
tion left is then how to find the explicit expression of the 
pseudo-rapidity e(0). This problem that has been recently 
solved by Caux and Konik [27]. We have to remind that 
on an interval L, each multi-particle state is in correspon- 
dence with a sequence of relative integers {N\, . . . , N n } 
(all different and here taken to be symmetric with respect 
to 0) entering the Bethe Ansatz equations [21, 22] 
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Thus any violation of the values of these ratios can be 
attributed to interaction, for the theory under study. 

Consider now a translation invariant initial state |^o) 
and as a post-quench Hamiltonian H the one of an in- 
teracting QIFT, with a two-body S'-matrix S(9). Ex- 
panding on the basis of the common eigenvectors of 
H and all higher charges, its general form is \ipo) = 



where ?(0) = — Hog S(0) is the phase-shift. Given the 
Ni's, the rapiditcs solution of these equations (the roots) 
are uniquely fixed and, for L — > oo, they form a distribu- 
tion p( r >(9). This function, in turns, determines through 
a linear integral equation the distribution p^ (9) of the 
holes, i.e. the solutions of eq.(10) corresponding to a 
skipped sequence of iVi's, since [21, 22] 
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J2p,n c P,n\Si,p---0 n ,p), with P p = m £,- smh i>p = 0. w h Grc ip(8) = d<,{9) / d9 . If one knew both p { - r \9) and 
Posing E p = mYd=i cosh^ )P , at any later time t the p (/l) (0), one could have access to e(0), since this quantity 



EV of a local observable 0(x, t) on \i))q) is 



is given by 
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(Z = (ipolipo)) ■ Taking the Dynamical Average, one ends 
up to the Diagonal Ensemble 
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As it stands, however, this expression is highly prob- 
lematic: all terms of the sum in the numerator as well 
as those present in Z are in fact divergent. Concerning 
those in Z, the divergencies come from the normalization 
of the eigenstates, (6' m , ...9[\9 1 ... 9 n ) = ]1. W - 9), 
which gives rise to [<5(0)] n when 9\ = 9. Concerning the 
numerator, the divergencies come from the Form Fac- 
tor (9' m . . . 0i|e>(O)||0i . . . 0„), once evaluated at 9[ = 9 t . 
These divergencies are an unavoidable consequence of the 
kincmatical pole structure of the Form Factors [28]. 

Luckily enough, the cure of this problem has been al- 
ready found by LeClair and Mussardo in the context of 
finite-temperature correlators [20] , where one encounters 
a similar difficulty: the remedy consists in employing the 



P M(0) ' 

In the thermodynamic limit there is a large number of 
quantum states compatible with the densities p^ and 
p^ h \ so there is a macroscopic entropy S[p^ r \ p^]. Usu- 
ally, the determination of p^ and p^ passes through 
the minimization procedure of the Free Energy [21, 22], 
a functional that in a QIFT involves all conserved charges 
[15, 26]. To fix the relative lagrangian multipliers may be 
not an easy task but, as shown in [27], one can actually 
bypass this problem and avoid altogether to follow the 
usual procedure adopting instead this plan of actions: (i) 
denoting by 9 P ^ the solution of eq.(10) relative to the 
generic mult i- particle state |0i ]P . . . 9 n ,p), associate to it 
the distribution p Pi „(0) = Yli $(@~@p,i) an d t a ke as den- 
sity of the roots of the initial state \i/jo} the weighted sum 
of all these distributions p^ (0) = J2 P \ c p,n\ 2 P P ,n{&) ', (ii) 
substitute this p^ into cq.(ll) and solve for p^ h \9) to 
get then e(0) from eq.(12). For the Generalized Free En- 
ergy one has F = -^ / d0 cosh log (1 + e _e W), while 
for the EV of the conserved charges on \ijjq) [26] 

<9e(0) 



Sto-i = ^ I d9 cosh 0/(0) 
2ir J 



da 2r , 




FIG. 3: Quench instances of {(: ip' ip :) versus the dimen- 
sionless coupling constant 7 of the Lieb-Liniger model. The 
thermal equilibrium curve is the dashed one. 



where 



/(*) = 



1 



e W + 1 



(13) 



is the filling fraction of the states, while dtjda^n-i is the 
solution of the linear integral equation 



de{6) 
9a 2n -i 



cosh[(2n-l)0]+ f ip(0-0')f(0') 



n de(0') 



da 2n -\ 



A nice check of eq.(3) is provided by the Lieb-Liniger 
model regarded as non-relativistic limit of the ShG model 
[23]. Using the Form Factors of this theory, one recov- 
ers the formulas obtained in [29] for the EV of the lo- 
cal composite fields : ip^ k (x)ip k (x) : of the Lieb-Liniger 
model, where the function f(p) employed in [29] is noth- 
ing else but the same f(0) of eq.(13). As an example, Fig. 
3 shows the EV of : tjv 2 (x)ip 2 (x) : versus 7, the dimen- 
sionless coupling of the Lieb-Liniger model, for different 
quench processes (this is the equivalent of Fig. 2 for an in- 
teracting case): no surprise, there is a spread of values of 
this quantity, since there are of course innumerable many 
ways to get out of equilibrium (and QIFT keeps memory 
of it) but only one way to be at equilibrium. 
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Note Added. There has been recently another proposal 
to compute the EV [30] , based on the Bethe Ansatz and 
checked for the free case of the quantum Ising model. 
Although very similar to mine, it is not clear to me that 
the two proposals are equivalent and further analysis is 
needed to prove or disprove their identity. 
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